IEEE Transactions on Biomedical Engineering

Non-invasive human-machine interfaces (HMIs) are critical in developing prosthetic systems that offer intuitive, simultaneous, and proportional control over multiple degrees of freedom (DOFs). This study introduces a novel system for intuitive concurrent control of hand and wrist movements using sonomyography based imaging of muscle activity. Our method uses a sparse set of ultrasound scanlines to reduce computational complexity while enhancing usability. We evaluated four regression techniques for wrist and hand angle prediction, focusing on performance with a reduced sonomyographic feature set. We also explored the feasibility of a sonomyography-based system by simulating various factors that could affect prediction, including feature selection and scanline count. Our findings demonstrate that Gaussian process regression excels in predicting wrist and hand angles with just eight equispaced transducers in offline settings. Real-time evaluations with 10 non-disabled participants showed a 93 % success rate for two-DOF tasks using linear regression. The system was tested with an individual with amputation, achieving a 46 % success rate for two-DOF control in a 2D space, even though the ground truth data for model training was collected from the contralateral limb. This study validates our sonomyography-based approach for accurate wrist and hand angle estimation, reducing complexity and demonstrating potential in real-world scenarios.

SignificanceDiffuse Correlation Spectroscopy (DCS) is an established technique for non-invasive monitoring of Cerebral Blood Flow (CBF), but existing applications primarily measure CBF changes averaged over cortical tissue volumes. A method capable of targeting blood flow within a particular intracranial artery would expand the utility of DCS for vessel-specific cerebral perfusion monitoring and continuous bedside assessment. AimWe aim to investigate the feasibility of targeted, non-invasive monitoring of blood flow in the Anterior Cerebral Artery (ACA) using DCS through optimization of probe geometry and placement. ApproachA custom-built DCS system operating at 785 nm was used to probe ACA from the glabellar region. Source-detector (SD) separation, probe orientation, and probe location were systematically optimized using lower-limb motor tasks and a mental arithmetic task. The optimized configuration was evaluated using an ACA-mimicking multilayer phantom and validated in forty healthy volunteers during lower-limb activation tasks and postural changes. ResultsAn SD separation of 17 mm, vertical probe orientation, and glabellar placement provided the highest sensitivity to ACA-related blood flow changes. Phantom experiments demonstrated sensitivity to flow changes in a vessel located 45 mm beneath the scalp and showed that the frontal sinus, cerebrospinal fluid, and the absence of cortical tissue beneath the glabella along the longitudinal fissure together provide the best optical window for probing deep ACA flow. Using the optimized probe configuration, significant increases in relative CBF were observed during standing leg marching (66.4 {+/-} 38.6%) and supine leg crunches (39.4 {+/-} 32.2%) (p < 0.01), with consistent responses during supine-to-stand postural transitions. ConclusionThe proposed DCS approach enables targeted measurement of blood flow within the ACA territory. This technique provides a framework for continuous, non-invasive monitoring of ACA perfusion and has potential applications in cerebrovascular monitoring and stroke care.

The features of the PPG (photoplethysmography) morphology are known to reflect age-related cardiac and vascular changes. In most contemporary wearables, PPG signals are acquired from distal sites such as the wrist and finger. The superficial temporal artery (STA), accessible at the temple region, is reached via a shorter arterial path from the aortic root than the radial circulation, and may therefore carry hemodynamic and aging information with less distance-dependent attenuation. We hypothesized that the morphology of the PPG at temple region (STA) would show stronger and more numerous age correlates than the PPG at the wrist. To test this, we extracted a common set of 89 pulse-morphology features, spanning raw-waveform timing/amplitude/area measures, ratios among them, derivative-based ratios, and spectral harmonic-ratio features. We compared an in-house temple-worn device which has PPG as one of the sensors, with a publicly available Microsoft Aurora-BP wrist-worn PPG dataset, and tested each feature's association with age. We identified 14 robust age correlates at the temple region, compared to 3 at the wrist. The temple's correlates spanned multiple morphological categories and showed a larger age-association than at the wrist. These results support the hypothesis that the temple region may be a more robust PPG measurement site than the wrist to extract age-related cardiovascular information, which motivates further investigation of temple-based cardiovascular sensing.

Muscle architecture shapes muscle function and changes with age, growth, training and disease, yet quantifying three-dimensional (3D) muscle architecture in vivo remains challenging. We introduce a hybrid fascicle tractography approach for freehand 3D ultrasound data that accurately reconstructs 3D muscle fascicles with respect to an objective, anatomically relevant coordinate system defined by the muscle's central aponeurosis. The hybrid approach combines Hessian-based fascicle detection with wavelet-based refinement to generate volumetric fascicle orientations. In a synthetic dataset with known ground truth, fascicle orientations and lengths were estimated with errors of [&le;]2{degrees} and ~1.5%, respectively. In vivo, the approach detected physiologically plausible fascicle lengthening in the human tibialis anterior following a passive plantar flexion rotation, whereas diffusion tensor imaging of the same muscle did not. The proposed method enables anatomically relevant, objective and non-invasive quantification of 3D muscle architecture in vivo, providing a practical framework for applications in clinical and applied muscle physiology.

Background and ObjectiveAccurately parameterizing dynamic, time-varying interactions in physiological systems is a methodological challenge, as global causal discovery methods may obscure transient, local fluctuations. This study introduces tempord, an open-source Python library designed to estimate local temporal orders and evaluate the short-term stability, directionality, and strength of causal links in non-stationary biological signals. MethodsThe algorithm estimates temporal relationships by keeping one signal stationary while iteratively shifting another one within a sliding window. To parameterize optimal inter-signal shifts (causal vector, CV), the framework utilizes linear modeling or time series distance metrics. The methodology was validated through a simulation study on synthetic bivariate signals with mathematically imposed dynamic phase delays, under both deterministic and noisy conditions. Furthermore, in-vivo capabilities were demonstrated by evaluating cardiorespiratory coupling dynamics across spontaneous and music-induced relaxation breathing states. ResultsThe simulation study demonstrated that the extracted CV trajectories precisely aligned with ground-truth temporal delays, assessed using mean absolute error and root mean square error for both noise-free and noisy synthetic data. In-vivo application demonstrated dynamic temporal stability and the detection of minor step changes during autonomic nervous system state transitions. ConclusionsThe tempord Python package bridges the gap between global causal discovery and local beat-by-beat statistical parameterization. It provides a robust "bottom-up" analytical instrument for investigating the transient mechanisms governing complex biological networks.

This paper presents a two-stage pipeline for implicit feature engineering in time series-based physiological stress detection using electrodermal activity (EDA) signals. In the first stage, we forecast three descriptive statistics of future EDA signals over short horizons (3, 5, and 10 seconds) based on a 60-second context window. In the second stage, a lightweight linear classifier detects stress from these predicted statistics. We evaluate three forecasting architectures spanning the domain expertise spectrum: a domain-specific bidirectional long short-term memory (BiLSTM) recurrent neural network, zero-shot and fine-tuned variants of Amazon Chronos T5 time series foundation model, and the Tabular Prior-data Fitted Network (TabPFN) applied to engineered physiological features. Experiments on the publicly available Wearable Stress and Affect Detection (WESAD) dataset, comprising chest-worn multimodal physiological signals from 15 subjects under baseline and stress conditions, demonstrate that the domain-specific BiLSTM achieves the highest classification performance, with area under the receiver operating characteristic curve (AUC) values ranging from 0.913 to 0.962. TabPFN follows with AUC values of 0.853-0.869, while Chronos variants yield 0.528-0.744. Notably, models using predicted features consistently outperform those using oracle features derived from the true future signals--the theoretical upper bound--suggesting effective noise filtering through learned sequence representations. Chronos models quickly reach performance saturation regardless of training depth, highlighting challenges in tokenizing continuous physiological time series. The proposed approach advances implicit feature engineering for wearable stress monitoring by leveraging forecasting as a powerful inductive bias, thereby improving robustness and providing insights into the limitations of the foundation model for physiological signals.

Background: Spinal cord injury (SCI) is frequently associated with orthostatic hypotension, defined by a sustained decrease in blood pressure upon assuming an upright posture due to impaired autonomic regulation. Cardiovascular spinal cord epidural stimulation (CV-scES) can regulate systolic blood pressure (SBP) in people with SCI, but stimulation paradigms are highly individualized. To make this treatment available to more patients, we developed an algorithm to tailor individualized CV-scES paradigms that closely mimic researcher-developed paradigms. Methods: We performed an offline analysis using datasets collected from eight individuals with SCI with epidural stimulators implanted over the lumbosacral spinal segments. During data collection, researchers modulated stimulation parameters with the goal of maintaining SBP between 110-120 mmHg. Each two-hour dataset included synchronized SBP and stimulation recordings. We ran optimization analyses offline to determine temporal requirements before modifying stimulation amplitude to mitigate out-of-range SBP. Results: The algorithm parameters that best matched researcher-selected stimulation changed relatively quickly during the first 12 min (one every ~40 sec), and more slowly thereafter (one every ~79 sec). Overall, algorithmic stimulation closely tracked researcher-controlled stimulation, with a mean correlation coefficient of 0.94. To evaluate online performance, we tested the algorithm in real time with a single participant. We found that a faster approach was needed to respond to changes in SBP caused by rapid, unpredictable events, such as postural changes. We implemented a sigmoid-based paradigm that determined the time to wait before changing stimulation as a function of the current SBP, with worse SBP values requiring faster responses. The new paradigm outperformed the original algorithm and researcher-controlled stimulation across measures of SBP stability, though recovery from a postural tilt maneuver remained slower than with researcher control. Conclusions: Our results indicate that algorithmic stimulation may minimize assistance required from researchers and participants, making CV-scES more feasible for clinical translation.

Transcranial low-intensity focused ultrasound (LIFU) offers a rare ability to noninvasively modulate human brain activity in superficial and deep neural structures. Yet the field still cannot reliably predict how a given acoustic protocol will affect a specific brain target. We synthesize 230 LIFU experiments targeting the human cerebrum from 136 articles reporting sufficient acoustic information and neuromodulatory outcomes. We compare associations between acoustic parameters and reported outcomes across cortical and subcortical targets. This reveals divergent parameter-effect landscapes across cortical and subcortical targets. Cortical outcomes follow a nonmonotonic organization across pulse repetition frequency and duty cycle, with distinct zones associated with facilitation, suppression, and null or ambiguous effects. Subcortical outcomes instead vary predominantly along duty cycle, shifting from suppression at lower values toward facilitation at higher values. Intensity and dose factors do not provide a consistent target-independent explanation for outcome polarity, and current evidence is insufficient to identify parameters promoting lasting effects. These findings challenge the assumption that acoustic protocols have generalizable effects across the brain and identify anatomical target as a central dimension that should inform LIFU protocol design. We provide SonoMap, an open, interactive web resource that supports a more systematic, evidence-based approach to human LIFU neuromodulation.

Surface electromyography (sEMG) offers a silent and wearable input modality, but its practical use is limited by variability across users and recording sessions. This study presents a compact CNN- Transformer model for decoding isolated handwritten digits from eight-channel sEMG signals. The model combines trainable signal preprocessing, convolutional feature extraction, and Transformerbased temporal modeling. It was evaluated on ten recordings from five participants using recordingseen classification, leave-one-recording-out (LORO) generalization, and few-shot adaptation. The model achieved a mean macro F1 score of 0.924 {+/-} 0.059 in the recording-seen setting and 0.619 {+/-} 0.252 under zero-shot LORO evaluation. Adaptation using two labeled trials per digit increased macro F1 to 0.828 {+/-} 0.112, while ten trials per digit achieved 0.925 {+/-} 0.053. The proposed architecture also outperformed classical and neural baselines in the controlled LORO benchmark. These results indicate that compact CNN-Transformer models, combined with lightweight target-recording calibration, provide a promising basis for adaptive sEMG-based input systems.

Marker-based optical motion capture remains the gold standard for kinematic analysis across biomechanical research. However, when tracking dynamic movements, severe occlusions often create missing data and noisy signals that require time-consuming manual cleaning. Existing automated methods usually rely on predefined skeleton models or need massive training datasets for machine learning solutions. This study presents the AMP: an automated marker based motion capture processing pipeline. The AMP is a fully analytical approach that uses hierarchical algorithms including 1) rigid-body constraint, 2) Kalman Filter state estimation and 3) kinematic pattern match to label markers and fill data gaps without manual human intervention. First, the algorithm identifies and removes mislabeling markers before any data scrubbing. During the data cleaning process, our hierarchical algorithm uses a translation-based step to track segments even when there are insufficient donor markers to form a full rigid body. To reconstruct missing marker trajectories, the pipeline initially uses Piecewise Cubic Hermite Interpolating Polynomial (PCHIP) interpolation for short marker gaps. However, the PCHIP creates massive spatial errors during long occlusions, the system automatically switches to rigid-body tracking for gaps over 40 ms. Results showed the pipeline drastically reduced processing time compared to manual human expert methods by 73.5%. In addition, the AMP kept spatial errors within 10 mm of the manual ground truth. By using this combined hierarchical approach to match marker processing accuracy, the pipeline provides a fast, general-purpose tool that fully automates data cleaning and saves hours of manual effort.

Carotid endarterectomy carries the risk of intraoperative cerebral ischemia, which is monitored by expert neurophysiologists through continuous electroencephalography (cEEG). Because expert availability is limited, we developed a hybrid novice-artificial intelligence (AI) system that detects ischemia using novice monitors with limited cEEG training. The hybrid system dynamically weights novice and AI inputs to arrive at a final output. Using four novices, we compared hybrid systems against experts alone, novices alone, and AI alone. Hybrid systems were statistically non-inferior to experts in sensitivity and false-positive rate (FPR), whereas novices alone were not. At 80% sensitivity, hybrid systems reduced FPR by half compared with the AI-only system, with similar benefits at 90% sensitivity. Further, the area under the precision-recall curve improved from 0.546 to 0.610-0.726, the area under the receiver operating characteristic curve improved from 0.957 to 0.967-0.971, and calibration improved compared with AI alone. These results highlight the potential of a hybrid system to monitor intraoperative cerebral ischemia.

Emerging MEMS-based audio devices generate ultrasonic acoustic output, but little is known about long-term biological effects of such exposures, which vary in frequency, intensity, duration, and coupling pathway. In this study, a novel transducer was characterized for its ultrasonic acoustic output, generating amplitude-modulated (AM) pressure pulses at a carrier frequency of approximately 200 kHz and operating frequency of approximately 100 kHz. To understand ultrasonic pressure wave transmission through the human middle and inner ear, sound pressure levels in the ear canal (PEC) and mechanical vibration velocities at the stapes (VST) and promontory (VPRM) were measured in human cadaveric temporal bones across frequencies up to 240 kHz. This represents the first such measurements at ultrasonic frequencies in human temporal bone specimens. Results demonstrated relatively consistent PEC measurements at ultrasonic frequencies across specimens, with minimal inter-specimen variation. At the operating frequency (~100 kHz), overall PEC was 82.2 dB SPL (SNR = 40.6 dB), while at the carrier frequency (~200 kHz), PEC increased to 90.6 dB SPL (SNR = 34.8 dB). VST and VPRM measurements at the operating frequency remained relatively independent of audio stimulus-drive frequency and voltage, with comparable magnitudes, suggesting that bone conduction pathways contribute significantly to inner ear ultrasonic exposure. When compared to existing safety guidelines, this device appears to meet specified criteria; however, given the limited scientific basis for these guidelines and unknown long-term effects, cautious application is recommended.

Accurate brain tumor segmentation from magnetic resonance imaging (MRI) remains a challenging task because supervised deep learning models require large quantities of annotated data, which are expensive and time-consuming to obtain. This study investigates whether synthetic MRI images generated using a Deep Convolutional Generative Adversarial Network (DCGAN) can improve U-Net-based brain tumor segmentation using synthetic data augmentation. Experiments were performed on the LGG-MRI dataset comprising 3,929 image-mask pairs. A baseline U-Net was first trained using the original training dataset. Synthetic MRI images were subsequently generated using a DCGAN, and threshold-derived pseudo-masks were assigned to the generated images to construct an augmented training dataset. The same U-Net architecture was then retrained using the augmented dataset and evaluated on an identical held-out test set. Compared with the baseline model, DCGAN-based augmentation increased the Dice coefficient from 0.2067 to 0.3037 and the Intersection over Union (IoU) from 0.1243 to 0.1918, while reducing the final test loss from 0.0474 to 0.0275. These results indicate that synthetic MRI augmentation was associated with improved segmentation performance under the experimental conditions of this study. However, the reliance on threshold-derived pseudo-labels and evaluation on a single dataset limit the generalizability of the results. The proposed workflow provides a reproducible implementation for evaluating DCGAN-based synthetic data augmentation in supervised brain tumor segmentation and establishes a baseline for future studies employing more reliable annotation strategies and broader experimental validation.

Neural temporal point processes (NTPPs) are powerful tools for modeling sequences of timestamped events with statistical temporal structure. Density-based NTPPs, in particular, are an interesting opportunity to merge the universal function approximation capability of neural networks with a defined statistical model in a way that has many potential applications. We demonstrate one such application to heartbeat dynamics, a physiologic point process. We specifically apply a lognormal mixture NTPP to compute instantaneous estimates of the mean and standard deviation of beat-to-beat intervals. We compare our results to the state of art (Barbieri et al.) point process model for heartbeat dynamics, which uses a more physiologically rigorous inverse Gaussian model. We find that the NTPP model maintains reasonable accuracy while improving upon robustness to noise.

Ultrafast power Doppler imaging combined with singular value decomposition (SVD) clutter filtering has become a standard approach for label-free microvascular ultrasound, enabling the visualization of small vessels without microbubble contrast agents. In the absence of contrast, however, SVD filtered Doppler images remain limited by the blur of the imaging system point spread function (PSF) and by a residual noise floor that reduces sensitivity at depth, which together hinder the resolution of fine microvasculature. Here we establish a sparse deconvolution framework to SVD filtered ultrafast power Doppler images. Each Doppler frame is processed in two cascaded stages: a Split-Bregman optimization that solves a regularized least-squares problem combining a sparsity prior and a Hessian continuity prior, followed by an accelerated Richardson Lucy deconvolution with an estimated system PSF. We first validated the framework on a simulation phantom with known ground truth, and then evaluated the framework on in vivo rat-brain plane-wave acquisitions obtained with a Verasonics Vantage system and a 15-MHz linear array. Compared to conventional SVD power Doppler, sparse deconvolution improved the resolution by around 4 and 8 times to lambda/2 and lambda/4, in the lateral and axial directions respectively. We further show that decomposing the Doppler signal into velocity bands before deconvolution disentangles slow and fast flow and yields a velocity-resolved microvascular map. Finally, applying the same framework to task-evoked functional ultrasound, we show that sparse deconvolution preserves the stimulus-locked cerebral-blood-volume response measured by conventional functional ultrasound while sharpening the corresponding activation map from a diffuse cortical region to discrete penetrating vessels. These results indicate that the sparsity-prior super-resolution principles established in label-free ultrafast Doppler ultrasound, and that sparse deconvolution can serve as a practical, contrast-free post-processing front-end for super-resolution microvascular and functional imaging.

Augmented reality glasses typically have poor optical efficiency (<1%), which is detrimental for high power applications, such as prosthetic vision. Maxwellian view optics are a candidate for increasing efficiency, but suffer from a very limited exit pupil and thus field of view. Many approaches to extend the eyebox inherently decrease the optical power. We present a design that uses a high-density VCSEL array and eye tracking to create an extended eyebox without sacrificing the advantages of a Maxwellian view display. Our system achieves a 20 degree field of view, a resolution 48.8 LP/mm (10 m), and an end-to-end optical efficiency of 83%

Conventional clinical assessment of upper limb kinematics relies on goniometry or visual estimation, methods that provide limited information on multi-segment 3D kinematics and can have poor inter-rater reliability. Inertial measurement units (IMUs) represent a viable alternative; however, extracting anatomically meaningful joint angles from skin-mounted sensors remains challenging owing to sensor-to-segment misalignment, soft-tissue artefact, and suboptimal orientation estimation under dynamic conditions. This study developed and validated a multi-segment IMU-based framework for upper limb kinematic assessment, addressing three stages of the measurement pipeline: (i) sensor fusion gain optimisation as a function of task mean segment angular velocity; (ii) accuracy of an anatomical landmark registration-based alignment method; and (iii) comparative evaluation of landmark-based alignment against manual and neutral-pose alignment approaches. Twenty participants with typical shoulder function performed clinical and functional upper limb tasks, recorded simultaneously with an optoelectronic motion capture system and six IMUs tracking the thorax, clavicle, scapula, humerus, forearm, and hand. Orientation estimation used a gradient-descent complementary filter with velocity-dependent gain modulation. Sensor-to-segment alignment was achieved via manual, neutral-pose, and landmark registration methods. Mean sensor orientation errors ranged from 1.5 to 3 degrees across segments/tasks and increased with angular velocity. Landmark registration produced substantially lower joint angle errors than manual and neutral-pose alignment across all joints (root mean square error range: 3.7-7.8 degrees; manual: 8.3-24.2 degrees; neutral: 15.5-35.5 degrees). These results indicate that IMU-derived kinematic accuracy is context-dependent, varying with joint and movement type, suggesting that accuracy should be established for the specific movements and joints relevant to each clinical application.

Vector flow mapping (VFM) reconstructs left-ventricular (LV) blood velocity from color-Doppler echocardiography by combining the measured beamwise component with physical and regularizing constraints. Analysis of the discrete VFM formulation shows that the inverse problem is intrinsically ill posed: the occurrence of singular modes can be predicted from the geometry of the segmented blood-pool domain, the imposed boundary conditions, and the degree of smoothing. These modes can propagate uncertainty along entire transverse bands of the reconstructed velocity field, yet conventional VFM neither quantifies this uncertainty nor allows for correcting the blood-pool segmentation. We introduce Bayesian VFM (B--VFM), a hierarchical framework that jointly infers radial and transverse velocities, a probabilistic blood-pool mask, their spatially resolved uncertainties, and hyperparameters weighting Doppler and segmentation fidelity, mass conservation, boundary conditions, and smoothness. The discretized posterior admits a closed-form gradient and exact Hessian, enabling computationally efficient, gradient-based MAP estimation, sampling, and direct analysis of ill-posed modes. Posterior inference combines Gibbs sampling of conjugate Gamma-distributed hyperparameters with conditional maximum-a-posteriori estimation and a Laplace approximation for the high-dimensional velocity and mask fields. To accommodate systematic departures from planar mass conservation, B-VFM can learn the covariance of the planar divergence residual from an ensemble of flows and incorporate it as a structured model-discrepancy prior. Independent chains converged reproducibly, while covariance priors learned from flow ensembles illustrated how model discrepancies can be incorporated into the inference. B--VFM was evaluated using Lamb-Chaplygin dipoles under ideal conditions and with Doppler corruption, Doppler voids, and segmentation defects, and using the Hicks-Moffatt family of spherical vortices to assess violations of planar mass conservation. The method produced smooth reconstructions, localized uncertainty near unreliable measurements and regions of model inconsistency, and used flow information to correct segmentation errors. Within the tested vortex family, the data-informed planar divergence prior reduced velocity bias and mask distortion. B--VFM thus provides an uncertainty-aware reconstruction method and a flexible foundation for future VFM formulations incorporating additional priors, observations, and physical models. Future work will evaluate the method using clinical data and more complex three-dimensional benchmark flows.

Cerebral hemodynamics are difficult to monitor continuously outside the laboratory. Optical head-worn wearables have been proposed for tracking cerebral blood-flow signals, but they require comparison with an established cerebrovascular reference before they can be interpreted. We evaluated a temple-worn optical wearable, Temple, that outputs a proprietary, dimensionless Brain Flow index, intended as a proxy for relative changes in cerebral hemodynamics, against transcranial Doppler (TCD) ultrasound, which measures blood-flow velocity in the middle cerebral artery (MCAv). Twenty-three healthy adults completed two physiological challenges that elicit distinct and acute cerebral hemodynamic responses: a cycle-ergometer exercise protocol and a stand-to-supine postural transition protocol. Twenty participants were analyzed per protocol. The Brain Flow index tracked MCAv in both protocols, with significant within-subject temporal correlations (median Pearson r = 0.795 and 0.799 for exercise and postural transition; p < 0.001) and directionally concordant, statistically significant transition responses for both increases and decreases in flow. Bland-Altman analysis of the normalized transition responses showed small mean biases between the two devices, consistent with similar relative response shapes. Because both signals were standardized within session before this comparison, it addresses the shape of the relative change rather than agreement in absolute units. The Brain Flow index reproduced the direction and time course of MCAv under both perturbations, including the postural transition, where heart rate moved in the opposite direction. Further studies using complementary modalities and additional cerebrovascular reactivity challenges are required to establish clinical use cases and cerebral specificity of the Brain Flow index.

ObjectiveElectroencephalography (EEG) is a key tool for studying auditory processing in cochlear implant (CI) users. In particular, EEG recordings obtained during continuous speech are becoming increasingly important for assessing speech and language processing in CI users, and may be utilized for neurofeedback. However, CIs also induce strong stimulation artifacts that are time-locked to the stimulus and mask the neural responses that have smaller magnitudes. Existing artifact reduction methods are typically based on event-related potentials (ERPs) or require manual component selection, making them unsuitable for naturalistic listening conditions or large datasets. ApproachWe develop CORSICA (CORrelation-baSed ICA artifact rejection), a reproducible, parameter-efficient method for CI artifact reduction in EEG responses to continuous speech. CORSICA operates on independent components (ICs) obtained through Infomax ICA and requires no manual component labelling, with performance governed by a single tunable threshold. It exploits the observation that CI artifacts temporally follow the audio signal without delay, whereas neural responses have an inherent lag due to auditory pathway latencies. For each IC, CORSICA computes the cross-correlation with the speech stimulus. Artifacts are identified by a high signal-to-noise ratio (SNR) of the correlation peak near zero lag, and the component is rejected if this SNR exceeds a threshold. To benchmark CORSICA, we evaluate two alternatives: a TRF-based SNR method, in which temporal response functions are fitted to each IC and artifact-driven peaks near zero lag are used for rejection, and a variant replacing ICA with second-order blind identification (SOBI) as the source separation step. Main resultsCORSICA effectively suppressed CI artifacts while preserving neural activity, enabling recovery of physiologically plausible TRFs with only 2% of ICs rejected. Both benchmark methods confirmed the validity of the SNR-based rejection framework, but CORSICA outperformed the TRF-based alternative in artifact suppression quality. Replacing ICA with SOBI as the source separation step required more ICs to be rejected, further supporting ICA as the preferred backbone for CORSICA. SignificanceCORSICA provides a fully objective, label-free approach to identifying CI artifacts in speech-evoked EEG data, with no manual intervention required. By centering artifact rejection on a single interpretable threshold, it offers a reproducible preprocessing standard for future EEG studies on speech processing in CI users. ConclusionOur findings demonstrate that objective CI artifact suppression in speech-evoked EEG data is feasible on the basis of the ICs temporal response patterns.